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Abstract. We derive the full phase diagram for a large family of exponential random graph 
models, each containing a first order transition curve ending in a critical point. 



1. Introduction 

We will treat a class of models of 'dense' random graphs, that is, simple graphs on n vertices in 
which the average number of edges is of order nP. More specifically we will consider exponential 
random graphs in which dependence between the random edges is defined through some finite 
graph, in imitation of the use of potential energy to provide dependence between particle 
states in a grand canonical ensemble of statistical physics. Exponential random graphs have 
been widely studied (see [U S] for a range of recent work) since the pioneering work on the 
independent case by Erdos and Renyi [2]. We will concentrate on the phenomenon of phase 
transitions which can emerge for dependent variables. Following analyses using mean-field 
and other uncontrolled approximations (see [HllO]) there has recently been important progress 
by Chatterjee and Diaconis [Ij, including the first rigorous proof of singular dependence on 
parameters. We will extend their result both in the class of models and parameter values under 
control and provide an appropriate formalism of phase structure for such models. 



2. Statement of Results 

We consider the class of models in which the probability of the simple graph G„ on n vertices 
is given by: 

where: Hi is an edge, H2 is any finite simple graph with p > 2 edges, ipn = i^niPi-ih) is the 
normalization constant, t{H^ Gn) is the density of graph homomorphisms H — > Gn'- 



|i/(G„)|ii'mi • 

and V{-) denotes the vertex set. Expectation of a real function of a random graph is denoted 
^01, main results are the following. 

Theorem 2.1. Restrict to > 0- Then the pointwise limit 

V'oo(/3i,/32) = hm M/^i,(^2) (3) 
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exists and is analytic in j3i and j32 off a certain curve /32 = which includes the end point 

/g^) = f 1 logfa - 1) - ^ , \ y (4) 



d d 

The derivatives CL^d -w^^oo o-f^ both discontinuous across the curve, except at (/3i,/32) 

opi op2 



where, however, all the second derivatives t—tttV'ooj tttttttt'^oo) o.nd 7— tttV'oo diverge. 

dj3f of3idf32 

Theorem 2.2. If in Theorem \2.1\ the graph H2 is a p-star, p > 2, the analogous result also 
holds for (32 < 0. 

Remark. A p-star has p edges meeting at a vertex. 

Corollary 2.3. The parameter space {(/3i, /32)} of each of the models consists of a single phase 
with a first order phase transition across the indicated curve and a second order phase transition 
at the critical point (/SJj/Sg)- 

To explain the language of phase transitions in Corollary 12.31 we first give a superficial in- 
troduction to the formalism of classical statistical mechanics within d-dimensional lattice gas 
models; for more details see for instance [8]. 

Assume each point in a d-dimensional cube 

C = {-n, -(n - 1), • • • 0, • • • (n - l),nY C Z'^ (5) 

is randomly occupied (by one particle) or not occupied, and assume there is a (many-body) 
potential energy of value a ^ associated with every occupied subset of C congruent to a 
certain H2 C C. The interaction is attractive if a < and repulsive if a > 0. We define the 
probability that the occupied sites in C are precisely c by 

-l3[,iEi{c)+aE2{c)] 

K'^ic) = ^ , (6) 

where the parameter /3 > is called the inverse temperature, the parameter G M is called the 
chemical potential, the normalization constant Z„(/3,^) is called the partition function. Hi and 
H2 are subsets of C with Hi a singleton and the cardinality \H2\ > 2, and Ej{c) is the number 
of copies of Hj in c. (We are using 'free' boundary conditions.) One of the basic features 
of the formalism is that the free energy density, -F„(/3,|u) = ln[Z„(/3, /i)]/n'^, contains all ways 
to interact with or influence the system, so that 'all' physically signiflcant quantities can be 
obtained by differentiating it with respect to j3 and //. For instance 

|-F.(/3,/x) = -/3E,,,{§}, (7) 

the (average) particle density. To model materials in thermal equilibrium, calculations in this 
formalism normally require that the system size be sufficiently large, and in practice one often 
resorts to using n — >■ cx). With this as motivation we tentatively define a 'phase' as a set of states 
(i.e. probability distributions) corresponding to a connected region of the (/3, /i) parameter 

space, which is maximal for the condition that lim — 

(/3,/i) are analytic in /3 and n 

n-j-oo dpWfi'^ 

for all j, k. Intuitively one associates a 'phase transition' with singularities which develop in 
some of these quantities as the system size diverges. An important simplification was proven by 
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Yang and Lee [H] who showed that the hmiting free energy density Foo(/3,/x) = hm 
always exists and that certain hmits commute: 

Qj+k Qj+k Qj+k 

This impUes that phases and phase transitions can be determined from the hmiting free energy 
density, and so a phase is commonly defined (see for instance [5]) as a connected region of the 
(/3,^) parameter space maximal for the condition that F^{(3,^) is analytic. 

Using the obvious analogues for random graphs, with /3i playing the role of — and (32 the 
role of —/3a (and therefore positive if and only if the model is 'attractive'), -0^ = tpniPi,/32) 
plays the key role of free energy density. We will show in Theorems 13.101 and 13.111 below that 
the limiting free energy density il^ooil^i-, h) exists, and the proof of Theorem 2 by Yang and Lee 
|llj . on the commutation of limits, then goes through without any difficulty in this setting, so 
we can again define phases and phase transitions through the limiting free energy density, as 
follows. 

Definition 2.4. A phase is a connected region, of the parameter space {(/3i,/32)}; maximal for 
the condition that the limiting free energy density, ■0oo(/3i, /32), is analytic. There is a j*^^ -order 
transition at a boundary point of a phase if at least one j''^ -order partial derivative of if) oo{Pi, 
is discontinuous there, while all lower order derivatives are continuous. 

Theorems I2.1l and l2.2l thereby justify our interpretation in Corollary 12.31 that each of our models 
consists of a single phase with a first order phase transition across the indicated curve except at 
the end (or 'critical') point (Pf,/]!^), where the transition is second order, superficially similar 
to the transition between liquid and gas in equilibrium materials. 

3. Proofs 

Chatterjee and Diaconis have proven that the main object of interest, ipooiPi, exists for 
all /?! and nonnegative /32 and is the solution of a certain optimization problem: 

Theorem 3.1 (Part of Theorem 4.1 in [l]). Fix one of our models and assume H2 has p > 2 
edges. Then for any —00 < /3i < 00 and < (^2 < 00, the pointwise limit V'oo(/?i) /?2) = 
lim '0n(/?i,/52) exists: 

n— ^00 

V'oo(A,/32)= sup f/?iu + ^2M^-^MlogM-^(l-n)log(l-n)). (9) 
?ie[o,i] V ^ ^ / 

The following detailed analysis of this maximization problem is therefore fundamental to 
understanding the phase structure of our models, and we now address it. 

Proposition 3.2. Fix an integer p > 2. Consider the maximization problem for 

/(n;/3i,/32) = /3in + /32nP-^nlogn-i(l-n)log(l-n) (10) 

on the interval [0, 1], where —00 < /3i < 00 and —00 < (^2 < 00 are parameters. Then there is a 
V-shaped region in the (/3i,/32) plane with corner point (/3f,^2) such that outside this region, l{u) 
has a unique local maximizer (hence global maximizer) u* ; whereas inside this region, l{u) always 
has exactly two local maximizers u\ and u^. Moreover, for every /3i inside this V-shaped region 
(/3i < f^l), there is a unique (52 = '7(/?i) such that the two local maximizers of l{u; j3i,q{f5i)) are 
both global maximizers. Furthermore q is a continuous and decreasing function of f3i . 
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Figure 1. The V-shaped region (with phase transition curve inside) in the 
plane. Graph drawn for p = 3. 

Remark. By the Lebesgue Differentiation Theorem, q being monotone guarantees that it is 
differentiable almost everywhere. 

Proof. The location of maximizers of l{u) on the interval [0, 1] are closely related to properties 
of its derivatives l'{u) and l"{u): 



r(n)=/3i+p/32uP-i-^log 



2 

1 

2u{l-u)' 

We first analyze properties of l"{u) on the interval [0, 1]. Consider instead the function 

1 



1"(U) = p{p - 1)P2UP-^ 

f l"iu] 
m[u) 



2p{p - l)nP-i(l - u) ' 



Simple optimization shows 



and the equality holds if and only if u = ^^-y-- Thus 



< {p-l)uP-^{l-u) < 
p-i 



p-l\P 



p^ 



< m{u) < oo. 



(11) 
(12) 

(13) 
(14) 

(15) 



2{p - 1)P 

The graph of m{u) is concave up with two ends both growing unbounded, and the global 



minimum is achieved at u = This implies that for (32 < w~-\\p^ ^"{'^) < on the whole 



p— 1 

interval [0,1]; whereas for P2 > 2(p-i)p ' '"(^) ^^^^ take on both positive and negative values, 



and we denote the transition points by ui and U2 {ui < ^^—^ < U2)- 
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Figure 2. Outside the V-shaped region, l{u) has a unique local maximizer 
(hence global maximizer) u*. Graph drawn for /3i = —0.8, /32 = 0.1, and p = 3. 




Based on properties of l"{u), we next analyze properties of l'{u) on the interval [0,1]. For 
1^2 < 2(p-i)p ' ^® monotonically decreasing. For f32 > 2(p-i)p ' ''(^) decreasing from to 
ui, increasing from ui to U2, then decreasing again from U2 to 1. 

Based on properties of l'{u) and l"{u), we analyze properties of l{u) on the interval [0,1]. 
Independent of the choice of parameters /3i and /32, l{u) is a bounded continuous function, 
l'{0) = oo, and /'(I) = — oo, so l{u) can not be maximized at or 1. For /32 < 2(p-i)p ' ^'(^) 
crosses the u-axis only once, going from positive to negative. Thus /(«) has a unique local 
maximizer ((hence global maximizer) u*. For /32 > 2(p-i)p i t'^s situation is more complicated 
and deserves a careful analysis. If l'{ui) > (resp. I'{u2) < 0), l{u) has a unique local maximizer 
(hence global maximizer) at a point u* > U2 (resp. u* < ui). If l'{ui) < < l'{u2), then l{u) 
has two local maximizers u* and with < ui < < n2 < 
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Figure 4. Along the lower bounding curve of the V-shaped region, I (u) has 
two zeros and U2, but only is the global maximizer for l{u). Graph drawn 
for Bi = -0.8, B2 = 0.769, and n = 3. 



U U 

Figure 5. Along the upper bounding curve of the V-shaped region, l'{u) has 
two zeros and U2, but only U2 is the global maximizer for l{u). Graph drawn 
for /3i = -0.8, 132 = 1.396, and p = 3. 



p— 1 

Notice that ui and U2 are solely determined by the choice of parameter /32 > 2{p~i)p ' ^^"^ 
vice versa. By p3|) . 

^'(^1) = ^1 + ^ - I (16) 

2(p — lj(l — til) 2 1 — ui 

1'{U2) = /3i + ^-^ r - ^ log (17) 

2(p - lj(l - if2) 2 I-U2 

1 1 ti 

n(n) = —, 71 r log . (18) 

^ ' 2(p- 1)(1 -u) 2 ^ l-u ^ ' 



Consider the function 
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It is not hard to see that n(0) = oo, n(l) = oo, n{u) is decreasing from to then increasing 
from to 1, and the global minimum value is 

"(^) = 25rb)-5'°^(^-')- 

This implies in particular that l'{ui) > for /3i > | log(p— 1) — 2(p-i) • '^^^ '-'^^^ possible region 
in the (/3i,/32) plane where l'{ui) < < /'(M2) is thus bounded by /3i < ^log{p — 1) — 2{p-i) 
and (32 > 2{p~i)p ■ 

We now analyze the behavior of l'{ui) and l'{u2) more closely when (3i and /32 are chosen 
from this region. Recall that by construction, ui < < U2- By monotonicity of n{u) on the 
intervals (0, and (^^, 1), there exist continuous functions a(/3i) and 6(/3i) of f3i, such that 
l'{ui) < for ui > a(/3i) and l'{u2) > for U2 > b{f3i). a(/3i) is an increasing function of f3i, 
whereas is a decreasing function, and they satisfy 

n(a(/3i))=n(6(/3i)) = -/3i. (20) 

Also, as /3i — )• —00, a(/3i) — ?• and — ?> 1. By p^ . the restrictions on and 7x2 yield 
restrictions on /32- We have l'{ui) < for /32 < m(a(/3i)) and Z'(n2) > for /32 > m(6(/3i)). 
Notice that m(a(/3i)) and m(6(/3i)) are both decreasing functions of and as /3i — )■ —00, they 
both grow unbounded. By construction, for every parameter value (/3i,/32), l'{u2) > l'{ui). 
Also, for fixed m(a(/3i)) is the value of /32 for which l'{ui) = 0, and m(6(/3i)) is the value 
for which l'{u2) = 0. Thus the curve m(6(/3i)) must lie below the curve m(a(/3i)). And 
together they generate the bounding curves of the V-shaped region in the (/3i,/32) plane where 
two local maximizers exist for l{u). It is not hard to see that the corner point is given by 

{(31, P^) = (i \og{p - 1) - . (See Figures 1-6.) 

Fixing an arbitrary (3i < P^, we examine the effect of varying P2 on the graph of l'{u). It 
is clear from (fTT|) that l'{u) shifts upward as (32 increases. As a result, as (32 gets large, the 
positive area bounded by the curve l'{u) increases, whereas the negative area decreases. By 
the fundamental theorem of calculus, the difference between the positive and negative areas is 
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the difference between /(U2) and l{ul), which goes from negative {l'{u2) = 0, is the global 
maximizer) to positive {l'{ui) = 0, U2 is the global maximizer) as /32 goes from m(6(/3i)) to 
m(a(/3i)). Thus there must be a unique /32 : m(6(/3i)) < /32 < m(a(/3i)) such that and U2 
are both global maximizers. We denote this /32 by q{(3i); see Figures [U and [H 

By analyzing the graph of l'{u), we see that the parameter values of are exactly 

the ones for which positive and negative areas bounded by l'{u) equal each other. An increase 
in /3i will induce an upward shift of l'{u), which must be balanced by a decrease in (^2 = qiPi)- 
Similarly, a decrease in /3i will induce a downward shift of l'{u), which must be balanced by an 
increase in f32 = q{/3i). This justifies that q is monotonically decreasing in Furthermore, 
the continuity of l'{u) as a function of /3i and /32 implies the continuity of (7 as a function of 
/3i. □ 

Corollary 3.3. The transition curve (32 = qWi) displays a universal asymptotic behavior as 
j3i -00: 

lim |g(/3i)+/3i| = 0. (21) 

Proof. By Proposition 13. 2^ it suffices to show that as f3i —00, /(u; /3i, — /3i) has two global 
maximizers and U2- This is easy when we realize that as /3i — )• —00, l{u; /3i, —/3i) — > —00 
for every u in (0,1). The limiting maximizers on [0,1] are thus = and U2 = 1, with 

l{ul) = 1{U*2) = 0. □ 

So far we have used results from [1] but have avoided specific reference to the framework of 
graph limits, developed by Lovasz et al, which was used to prove those results. We now need 
to refer directly to graph limits; for the notation and an introduction to this material see for 
instance [7j or [1]. 

Theorem 3.4. Let Gn be a random graph on n vertices in one of our models. For parameter 
values of (/3i,/32) in the upper half-plane (^2 > — ^^^^^rjy, the behavior of Gn in the large n limit 
is as follows: 

min (5n(G„, ii) — )• in probability as n 00, (22) 

ueU 



where U is the set of maximizers of ilO\) . 

Proof. The assumptions of Theorems 4.2 and 6.1 in [T] are satisfied for parameter values /32 > 
~ p(p^i) ■ By Proposition 13.21 along the curve (/3i, g(/?i)), the maximization problem (fT0]l is 
solved at two values u* and U2; whereas off this curve, it is solved at a unique value u* . Thus 
in the large n limit, along the curve (/3i, g(/3i)), G„ behaves like an Erdos-Renyi graph G{n,u) 
(u picked randomly from u* and Ug); whereas off this curve, G„ is indistinguishable from the 
Erdos-Renyi graph G(n, M*). □ 

Corollary 3.5. Fix any (32 > (^2- ^ edge, so t{H,Gn) is the edge density of Gn- 

Then there exists a continuous and decreasing function q~^{(^2) such that 

lim F^y^t{H,Gn) >U2) = 1 if Pi > q-H(32) (23) 

n— >oo 

and 

lim F^^'^'{t{H,Gn) < ni) = 1 if (3^ < q'^h)- (24) 
Here ui and U2 are defined as in the proof of Proposition \3. 2[ m{ui) = m{u2) = (^2- 
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Remark. ^45 /32 — > oo, ui ^ and U2 ^ 1 and the jump is noticeable even for relatively small 
values of 1^2 ■ 

Proof. As g(/3i) is a continuous and decreasing function of the inverse function g~^(/32) 
exists and is also continuous and decreasing. We examine the effect of varying /3i on the 
graph of l'{u) (and hence on the global maximizers of l{u)). First note that varying j3i does 
not change the shape of l'{u). Inside the V-shaped region, there are three cases. Recall that 
* < ui < 2^ < U2 < For /?! = (7~^(/32), positive and negative areas bounded by l'{u) 



equal each other, thus and U2 are both global maximizers. For (3i < q~^{f32), the graph 
of l'{u) shifts downward, negative area exceeds positive area, thus is the global maximizer. 
For Pi > q~^{/32), the graph of l'{u) shifts upward, positive area exceeds negative area, thus 
U2 is the global maximizer. Outside the V-shaped region, there are two cases. Below the lower 
bounding curve, l'{u) has a unique local maximizer u* < ui. Above the upper bounding curve, 
l'{u) has a unique local maximizer u* > U2. Our conclusion then follows from Theorem 13.41 □ 

Theorem 3.6. Assume that in one of our models H2 is a p-star (p > 2). For all parameter 
values (/3i,/32), the behavior of Gn in the large n limit is as follows: 

min (5n(G'„, tt) — > in probability as n —?> 00, (25) 
where U is the set of maximizers of lilU\) . 

Proof. This follows from related results in [1]. We separate the parameter plane {(/3i,/32)} 
into upper and lower half-planes. The upper half-plane (/32 > 0) satisfies the assumptions of 
Theorem 4.2, and the lower half-plane (/32 < 0) satisfies the assumptions of Theorem 6.4. By 
similar reasoning as in Theorem 13.41 the rest of the proof follows. □ 

Real and complex analyticity are both defined in terms of convergent power series. Any 
problem in the real analytic category may be complexified and thereby turned into a complex 
analytic one, and any complex analytic situation with real coefficients is obviously real analytic 
and can thus be treated with real analytic techniques. The following analytic implicit function 
theorem may be interpreted in either the real or complex setting. 

Theorem 3.7 (Krantz-Parks [6]). Suppose that the power series 

F(x,y) = ^a,,fcxV (26) 

a,k 

is absolutely convergent for \x\ < Ri, \y\ < R2. //ao,o = and oq^i ^ 0, then there exist rp > 
and a power series 

fix) = cax- (27) 

|a|>0 

such that (21) is absolutely convergent for \x\ < tq and F{x,f{x)) = 0. 



Proposition 3.8. Off the end point il3i,f32), the local maximizer u* for l{u; /3i, /32) (ui and u\ 
if inside the V-shaped region) is an analytic function of the parameters f3i and (32- 

Proof. A local maximizer u* for l{u) is a zero for l'{u) with the additional property that l'{u) 
would change sign from positive to negative across u = u* . Fix a choice of parameters (/3J, /32) 7^ 
(/3J,/3^). Set X = (/3i - /3[,/32 - /3^) and y = u - n*(/3;,/3^). The function I' {u; f3i, f32) is thus 
transformed into a function F{x,y). It is clear that /'(n; /32) is analytic for u G (0,1), 
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/?! G (—00,00), and (^2 £ (—00,00). Recall that l{u) bounded continuous, /'(O) = 00, and 
/'(I) = —00 implies that u* can not be or 1. It follows that the transformed function 
F[x^ y) has the desired domain of analyticity, and is locally absolutely convergent. As for the 
coefficients, 

ao,o = F(0,0) = r(n*(/3l,/3^);/3;,/3^) = (28) 

by construction, and 

dF 

ao,i = ^(0,0) = l"iu*W[,P'2);P[,P2) / 0, (29) 

as can be seen from the proof of Proposition 13. 2i In more detail, for < f^!^, l"{u) is positive 
for all u in (0, 1). And for = l"{u) = only for u = which coincides with u* only 
for P'^ = 13^. Lastly, for > /32, l"{u) = only for u = ui and u = U2- Three possible 
situations might occur. Outside the V-shaped region, l{u) has a unique local maximizer u*, 
with u* > U2 if above or along the upper bounding curve, and u* < ui if below or along 
the lower bounding curve. Inside the V-shaped region, l{u) has two local maximizers u'^ and 
U2, with ti* < ui < 2^ < U2 < Ug. All the conditions of Theorem 13.71 are satisfied, thus 
f{x) = u*(fii,P2) - w*(/3l,/32) converges for (/3i,/32) close to {/3[,(3'2). □ 

Proposition 3.9. Off the phase transition curve, l{u*) = max/(n; /32) (Kul) or l{u2) if 
inside the V-shaped region) is an analytic function of the parameters (3i and (32. 

Proof. It is clear that Z(m;/3i,/32) is analytic for u € (0,1), /3i G (—00,00), and /32 G (—00,00). 
Outside the V-shaped region, l{u) has a unique local maximizer u* in (0, 1), which is analytic 
in /?! and (32 by Proposition 13.81 Inside the V-shaped region, l{u) has two local maximizers li* 
and U2, both have values in (0, 1) and are analytic in (3i and (32 by Proposition 13.81 Below the 
phase transition curve, maxZ(u) is given by l{ul), which coincides with l{u*) along the lower 
bounding curve. Above the phase transition curve, max/(ii) is given by l{u2), which coincides 
with l{u*) along the upper bounding curve. Our claim follows by realizing that compositions 
of analytic functions are analytic as long as the domains and ranges match up. □ 

Theorem 3.10. Let Gn be a random graph on n vertices in one of our models. The limiting 
free energy density ip^o = lim„_>.oo is an analytic function of the parameters (3i and (32 off 
the phase transition curve in the upper half-plane (32 > ~ p(p^i) ■ 

Proof. The assumptions of Theorems 4.2 and 6.1 in [Ij are satisfied for parameter values (32 > 
~ p(J'-i) ■ O^^^ claim then follows from Proposition 13.91 □ 

Theorem 3.11. Assume that in one of our models H2 is a p-star (p > 2). The limiting free 
energy density ip^o = hnin^oo V'n o,n analytic function of the parameters (3i and (32 off the 
phase transition curve. 

Proof. The assumptions of Theorems 4.2 and 6.4 in [T] are satisfied. Our claim again follows 
from Proposition 13.91 □ 

Lemma 3.12 (Lovasz-Szegedy [7]). Let U,W : [0,1]^ — ?> [0,1] be two symmetric integrahle 
functions. Then for every finite simple graph F, 

\t{F,U)-t{F,W)\ < \EiF)\-5n{U,W). (30) 



Proof of Theorems \2.1\ and \2.SX The stated analyticity is proven in Theorems 13.101 and 13.111 
so we only need to examine the situation along the phase transition curve. We know from 
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Theorems 13.41 and 13.61 that G„ converges in probabihty to u*, off the curve. By Lemma 13.121 
t{Hi,Gn) then converges in probabihty to t{Hi,u*). As t{Hi,Gn) is uniformly bounded in n, 
this imphes that 

Ef^^^P2{\t{Hi,Gn) -t{Hi,u*)\} ^ as n ^ oo. (31) 

Therefore 

E0,,lsMHuGn)} ^ E^,,f^MHi,u*)} = n*(/3i,/32) = t^V'oo(/3i, /32) as n ^ oo. (32) 

opi 

Similarly, 

E0,,fs,{tiH2,Gn)} ^ Ep,,p,{t{H2,u*)} = {u*{f3i,l32)r = -^i^UPi. P2) as n ^ 00. (33) 

d d 

By Corollary 13.51 these two first derivatives -rrpri'oo and ttttV^oo are discontinuous across the 

curve (except at the end point). Let us now take a closer look at the behavior of -000 at the 
critical point. Recall that l'{u] f3i,f32) is monotonically decreasing on [0, 1], and the unique zero 
is achieved at 2^. Take any < e < i. Set S = min{/'(2^ - e),-l'{E^ + e)}. Consider 
{(^17^2) so close to (/Jfj/Sg) such that \j3i — /3J| +p\(32 — /^^l < 6. For every u in [0, 1], we then 
have \l'{u; Pi, 132) — l'{u] fSf, /Sf)! < 5. It follows that the zeros n*(/3i, /32) (n| and U2 if inside the 

d 

< e, which easily implies the continuity of ttt-V'oo 
d 

and -^r-z-ipoo at (/Sf', ^2)- To see that the transition at the critical point is second-order, we check 

OP2 

the second derivatives of V'oo in its neighborhood. Off the phase transition curve, 

d 1 
n-s.00 dpf opf opi l"{u*) 



V-shaped region) must satisfy 



p-i 
p 



52 , 92 , a , (p(n*)P-i)2 



But as was explained in Proposition [321 l"{u*) converges to zero as (/3i, /32) approaches (/3J, /S^); 
the desired singularity is thus justified. □ 

4. Summary. 

Much of the literature on phase transitions in exponential random graph models uses tech- 
niques, such as mean-field approximations, which are mathematically uncontrolled. As such 
they have been useful in discovering interesting behavior but they can be misleading in detail. 
For instance, although phase transitions have been discovered in this way for the two-star {H2 a 
two-star) p] and edge-triangle {H2 a triangle) models [lOj, the approximation leads to an error 
in the qualitative nature of the transition, attributing phase coexistence to the full V-shaped 
region of Figure [1] rather than just the curve P2 = qWi)] in other words it does not distinguish 
the local maxima in the region from the global maxima. 

Chatterjee and Diaconis flj gave the first rigorous proof of singular behavior in an exponen- 
tial random graph model, the edge-triangle model. Our paper is an extension of this important 
first step; besides extending the models and parameters under control we have provided a math- 
ematical framework of 'phases' which we hope will be useful in motivating future mathematical 
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work in this subject. Our results show that all models with 'attraction' {P2 > 0), and also p-star 
models with 'repulsion' {P2 < 0) exhibit a transition qualitatively like the gas/liquid transition: 
a first order transition corresponding to a discontinuity in density, with a second order criticial 
point. In Theorem 7.1 of [l] Chatterjee and Diaconis suggest that, quite generally, models with 
repulsion exhibit a transition qualitatively like the solid/fluid transition, in which one phase 
has nontrivial structure, as distinguished from the 'disordered' Erdos-Renyi graphs, which have 
independent edges. We have not yet been able to extend our results to this regime, which is an 
important open problem. 
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